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ABSTRACT 

We present near-infrared H2, radio CO, and thermal infrared observations of the nearby massive star-forming 
region Cepheus A (Cep A). From H2 bow shocks arranged along four distinct jet axes, we infer that the massive 
protostellar source HW2 drives a pulsed, precessing jet that has changed its orientation by about 45° in roughly 
10'* years. The current HW2 radio jet represents the most recent event in this time series of eruptions. This sce- 
nario is consistent with the recent discovery of a disk around HW2, perpendicular to the current jet orientation, 
and with the presence of companions at projected distances comparable to the disk radius. We propose that the 
Cep A system formed by the disk-assisted capture of a sibling star by HW2. We present a numerical model of a 
15 M© star with a circumstellar disk, orbited by a companion in an inclined, eccentric orbit. Close passages of 
the companion through or near the disk result in periods of enhanced accretion and mass loss, as well as forced 
precession of the disk and associated orientation changes in the jet. The observations reveal a second powerful 
outflow that emerges from radio source HW3c or HW3d. This flow is associated with blueshifted CO emission 
and a faint H2 bow shock to the east, and with HH 168 to the west. A collision between the flows from HW2 
and HW3c/d may be responsible for X-ray and radio continuum emission in Cep A West. 
Subject headings: ISM: jets and outflows — ISM: Herbig-Haro objects — ISM: kinematics and dynamics — 
ISM: individual: CepA — stars: formation 



1. INTRODUCTION 

Collimated, bipolar outflows accompany the birth of young 
stars from the earliest sta ges of star formation to the end of 
their accretion phase (e.g. lReipurth & Ballyll200lllBally et alJ 
I2OO6I) . The structure and kinematics of these flows provide a 
fossil record of the mass-loss histories of the associated young 
stellar objects (YSOs); the most distant shocks trace the old- 
est major mass-loss events while inner jets and shocks provide 
clues about the recent processes. The structure and symme- 
tries of outflows record orientation changes of the underlying 
accretion disk and motion of the outflow sources relative to 
the surrounding interstellar medium (ISM). For example, C- 
shaped outflows may indicate deflection by a side -wind (e.g., 
within expanding plasma of the Orion Nebula, iBallv et all 
I2OO6I) : other C-shaped bends may result from motion of the 
source YSO through the medium, possibly ejected by dynam- 
ical processes in a youn g cluster (e.g., HH 49 8 in NGC 1333 , 
|BalliL& Reipurth''2001'; HH 366 in Barnard 5, Yu et al. 199Sj). 
Point-symmetric S- or Z-shaped bends can be produced by the 
precession of a circumstellar disk. Examples of this kind of 
symmetry include HH 199 in the L12 28 cloud dBallv et alj 
Il995h and 20126H-4104 (ISu et al.ll200'7l) . Such outflow orien- 
tation changes may provide indirect evidence for a companion 
star in an orbit inclined relative to the disk, inducing forced 
precession. 

In this paper, we present images of shocked outflows in the 
Cep A star-forming region, together with radial velocity maps 
of CO emission. We interpret the outflow morphology as ev- 
idence that the massive protostar HW2 drives a pulsed, pre- 
cessing jet whose orientation changes may be induced by the 
periastron passage of a moderate-mass companion in an ec- 



centric, non-coplanar orbit. We present a numerical model 
that demonstrates the plausibility of jet precession forced by 
a companion for the case of HW2. 

1.1. Overview of Cep A 

The Cep A star-forming complex contains the second near- 
est region of massive star formation, after the Orion com- 
plex. Located at a distan ce of 725 pc ( B laauw et al.l 1 1 9591 : 
ICrawford & Barnesll970h , the Cepheus OB3 association con- 
tains a 20 by 60 pc molecular cloud that houses six loc alized 
peaks of CO emission, designated Cep A through F (SargenJ 
1 197 71 |1979(). Cep A contains dense molecular cl umps 
J^rrelles et al. 1993 ), mole cular outflows (R odrigue z et alJ 
fl980a:,Naravanan & mikeiill996l:lG6mez et aL,1999ft . H2O 



and OH masers jCohen et al.l ] 1 984 ), hyper-compact H II 
regions (iHughes & Wouterloot 1984), variable radio con- 
tinuum sources (Garayetal. 1996), Herbig-Haro objects 
(Hartigan et al. 1986), bright shock-excited H2 emission 
(iHartigan et aL,2000) , a cluster of far-infrared (FIR) sources 
with a luminosity of 2.5 x 10"* Lq (Koppenaal et al. 1979; 



Electronic address: ncunningham2@unl.edu 

Electronic address: Nickolas. Moeckel (3* colorado.edu I 



Electronic address: John.Bally@colorado.edir| ~ 



and a cluster of Class I and Class II YSOs ( Gutermuth et 
l2005h . The bulk of the region's lu minosity Ukely arises 
from radio sources HW2 and HW3c/d (iHughes & Wouterloot! 
[T984b that are associated with bright H2O masers. 

Cep A contains a massive b ipolar molecula r outflow 
aligned primarily east-west dRodrfguez et al.l Il980bl) , 
but with additional cor npone nts aligned no rt heast-- 
southwest ( Ballv & Land 119901: iTorrelles et alj 119931: 
[Narayanan & Walker.. 1 9961 iGomez et al.ll 19991) . Tlie central 
2' region contains high velocity (HV) as well as more 
compact extremely high velocity (EHV) CO components 
with radial velocities ranging from -50 to 70 km s"' relative 
to the CO centroid (Nara yanan & Walkeilll99 6). The axis of 
the EHV outflow is rotated roughly 40° clockwise relative 
to the HV outflow on the plane of the sky. The smaller 
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spatial extent together with the higher velocity suggests that 
the EHV flow traces a younger outflow component. Several 
self-absorption dips follow trends seen in the low-velocity 
line wings, with regions east of HW2 blueshifted and west of 
HW2 redshifted. Thus cooler, self-absorbing gas traces the 
low-velocity bipolar outflow rather than the quasi-stationary 
ambient medium. At low velocities, there are additional 
blue- and redshifted components centered on HW2 that are 
oriented northeast-southwest. 

The Cep A outflow complex contains several Herbig-Haro 
objects, including the extremely bright HH 168 located about 
90" due west of HW2, and several fa inter bow shocks lo- 
cated to the east (iHartigan et alJl2000l) . Fainter HH objects 
(HH 169 and 174) are located in the eastern, blueshifted lobe. 
Near-inrared (NIR) images show an extremely bright reflec- 
tion nebula centered on HW2 with an illumination cone that 
opens toward the northeast. The 2.12 fim H2 line exhibits a 
complex, filamentary structure. 

There are multiple luminous sources in the complex of HW 
radio sources at the core of the Cep A region, led by HW2, 
and including HW3c, HW3d, and perhaps HW3a, HW8, and 
HW9. To the west, it has been proposed that radio source W-2 
may also be internally heated. Because some of these (HW2, 
HW3d, W-2) have the radio signatures of jets, these are likely 
candidates for driving the outflows observed in CO and other 
molecules as well as the shocks traced by H2, Fe II, and optical 
HH emission. The complexity of the outflows, their multiple 
orientations, ages, and velocities, and the unclear morphology 
of some of the shock features make Cep A challe nging to in- 
terpret. Are deflections involved, as suggested by Goetz et al.l 
([1998)? Are there additional outflow sources we have not yet 
detected? Are we confusing externally shocked or jet-heated 
features with self-luminous sources? 
HW2 is the strongest radio continuum source in the re- 
ion wifli a flux d ensity of 15.8 ± 0.3 mJy at 14.9 GHz 



Table 1 

NIC-FPS observations of Cepheus A 



iGaray et anil996l) : it may be the most luminous source in 
Cep A with L « 10"^ Lq, implying a mass of 15-20 Mq. 
HW2 is located near the center of a small (< 5") , dense 
clump apparent in the NH3 maps of iTorrelles et al.l (1 1993b 
at the tip of one of the larger NH3 concentrations. Multi- 
frequency observations by Rodriguez et al. (1994) using the 
VLA show that HW2 is elongated, and that both its spectral 
index and size versus frequency relation match those expected 
for a bi conical thermal jet. Higher angular resolution obser- 
vations ('Hughe s et al.lll995l:lHoare & Ga rrington'1995) show 
that HW2 contains a clumpy radio continuum jet whose knots 
have proper motions of order 500 km s"' (ICuriel et al.li2006l) . 
Many models of the region postulate that HW2 is the major 
or sole source of outflows in the region. 

The recent interferometric observations of the dust con- 
tinuum, free-free emission, and several molecular tracers 
with arcsecond and subarcsecond angular resolution have 
shown that HW2 is surrounded by a complex dust distribu- 
tion and several close-by protostars. A hot core that must 
contain a mod erate-mass YSO is loc ated about 400 AU east 
of HW2 (Mart m-Pintado et al.ll2005l) . VLA observations re- 
veal a supposedly low-mass protostar that emits radio waves 
an d is located at t he cent er of the H2O maser arc detected 
by ITorrelles et all (l2001bh . HW2 is surrounded by a cir- 
cumstellar disk at least several hundred AU in radius and 
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2005 Dec 1 1 


H2-2.12ftm 


10 


240 


2400 




Ks 
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20 


100 


2006 Jan 12 


H2-2.12^tm 
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240 


1200 




Ks 


20 


12 


240 


2007 Jan 29 
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15 


20 


300 




H 


15 


20 


300 




Ks 


15 


20 


300 



likely by s e veral l ow-mass proto stars ( Rodriguez e t al. 2005; 
Pateletal. 2005; Broganetal. 2007; Co mito et all I2007t 
Jimenez-Serra e t al. 2007; Torrelles et al, 2OO7I) ! 
The radio sources Cep A HW3d and HW3b are also associ- 



ated with a clusters of H2O masers. HW3c and/or HW3d are 
suspected to be moderate-mass protostars. They are located 
5" south of HW2 and close to the western end of a nearly 
continuous chain of radio sources along the southern rim of 
the Cep A East radio source complex. The eastern part of this 
chain may in part trace a radio continuum jet from HW3c/d at 
about RA. « 100°. The counter-jet direction points directly 
at the radio source W-2 at the eastern end of HH 168 in Cep A 
West. 

Su bmillimeter wavelen gth continuum interferometry at 875 
/zm ( iBrogan et alj l2007h shows strong dust emission from 
HW2, HW3c, and extended emission to the southwest near 
HW3a. No submm continuum was detected from HW3d, 
leading these authors to conclude that HW3c is most likely 
to harbor the second most luminous and massive YSO in the 
Cep A core. In this interpretation, radio source 3d may trace 
part of a thermal jet from HW3c. The associated maser emis- 
sion may be an indicator of shocks in a dense molecular gas. 
For the rest of this paper, we will assume that the dominant 
energy source is embedded in HW3c, but our interpretation 
does not depend critically on which peak, HW3c or 3d, con- 
tains the source. 

The extinction to the HW sources is extremely large with 
Ay « 500 to 1000 magnitudes. Thus, none of these radio 
continuum sources is visible in the NIR. An extremely bright 
infrared (IR) reflection nebula emerges from the vicinity of 
HW2 toward RA. w 60°. The IR continuum polarization vec- 
tors indicate that the illumin ation is coming from the d irection 
of HW2, HW3, and HW8 dColome & HarvevI [19951) . How- 
ever, the accuracy of polarization measurements is insufficient 
to determine the source unequivocally. 

2. OBSERVATIONS 
2.1. NIR Imaging 

We observed Cep A on three nights from 2005 to 2007, 
using the Near Infrared Camera/Fabry Perot Spectrograph 
(NIC-FPS) on the Astrophysical Research Consortium 3.5m 
telescope at Apache Point Observatory, New Mexico. The 
NIR emission in the Cep A region spans 9' east to west, and 
requires at least two telescope pointings for complete cover- 
age with the 4.6' x4. 6' field-of-view of NIC-FPS. On each 
occasion, the total exposure time was evenly split between a 
western field and an eastern field. We obtained broadband im- 
ages in J, H, and filters of the Mauna Kea filter set, and in 
a narrowband (0.4% bandpass) filter centered on the H2 2.12 
/im line. The dates and details of observations taken on each 
night are listed in Table[T] 

Images were taken in dithered sets of 5 or 15 frames; indi- 
vidual frames were reduced by performing dark and sky sub- 
traction, flat fielding, and correction for instrumental geomet- 
ric distortion. Frames in each set were then co-aligned and 
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Figure 1. The 2.12 /im H2 emission (green) superimposed on emission from high-velocity CO with —27 < Vlsr < km s ' (blue) and —3 < Vlsr < 6 km s ' 
(red). In this (and the next two figures) the purple circle shows a beam-sized spot at the location of HW2. 
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Figure 2. The 2.12 /im H2 emission (green), and emission from moderate-velocity CO with -18 < Vlsr < ~13 km s (blue) and —6 < VLsR < ~3 km s (red). 
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median combined to eliminate bad pixels and to reduce noise. 
Registration to absolute coordinates was achieved by match- 
ing locations of many point sources in each combined image 
to stellar J2000 coordinates for the region from the 2MASS 
project. Once properly registered to this common frame, all 
the images taken in a single filter were co-added; each region 
of the combined image was scaled to properly account for a 
varying number of overlapping images and differing exposure 
times. 

As the H2 2.12 fim narrowband filter falls within the band- 
pass of the Ks filter, the latter was used as a proxy for the con- 
tinuum that contaminates the H2 emission in the narrowband 
image. We removed the continuum component by subtracting 
a scaled version of the final Ks image from the same nights 
as the H2 observations. The scale factor was chosen based 
on the best continuum removal as judged by eye in the sub- 
tracted image. Variations in seeing between the narrowband 
and broadband exposures result in the incomplete removal of 
some stars. 

2.2. Optical spectroscopy 

Visual wavelength long-slit spectra of the shocks located 
east and northeast of Cep A HW2 were obtained with the 
Double Imaging Spectrograph (DIS) at the f/10 Nasmyth fo- 
cus of the 3.5 meter reflector at the Apache Point Observatory 
on 2007 January 9. A 1.5" wide by 5' long slit was used with 
a grating providing a spectral resolution R = A/AA w 5,000. 
At each slit orientation three images with exposure times of 
10 minutes were obtained. 

2.3. Radio observations 

The Nobeyama Radio Observatory (NRO) 45 m radio tele- 
scope was used in February 1987 to map the distribution of 
J=l-0 '^CO from the Cep A outflow complex with 15" res- 
olution at 115 GHz. Nobeyama spectra were collected on a 
uniform grid with a 10" spacing over a 5 by 9 arcmin re- 
gion containing the Cep A outflow complex. Over 1900 in- 
dividual positions were observed for an average integration 
time of about 1 min per position, with 20 seconds on source 
and 20 seconds off source plus 10 seconds of dead time used 
for reading the spectrometer, and moving the antenna after 
each 20 second integration period. Approximately 14 sepa- 
rate 10-hr periods were devoted to observing Cep A. A cooled 
Schottky receiver with a single sideband receiver temperature 
TssB ~ 350K was used in conjunction with both a high (25 
kHz/channel) and a low resolution (100 kHz/channel) 2048 
channel acousto-optic spectrometer (AOS). All data were cal- 
ibrated using the standard chopper wheel method. Pointing 
was monitored every 2 hr by observations of 43 GHz SiO 
maser emission from nearby late-type giant stars using a sep- 
arate receiver system. Both high- and low-resolution AOSs 
were used to record spectra. The noise temperature of this re- 
ceiver fluctuated between Tssb ~ 800K and infinity due to a 
problem, so data in these lines are available only over parts of 
the source. Spectra were obtained by position-switching to an 
emission-free reference position about 20' to the northwest. 
The data were reduced using the COMB package developed 
by R. W. Wilson at AT&T Bell Laboratories. 

2.4. Thermal IR Imaging 

Observations of the Cep A region were made using the 
Keck Observatory facility mid-IR camera Long Wave Spec- 
trometer (LWS) on 2002 November 16 November LWS is 



a mid-IR imaging and spectroscopy instrument mounted on 
the forward Cassegrain focus of Keck I, employing a Boe- 
ing 1 28 X 128 As: Si BIB ar ray with a 10.2" x 10.2" field of 
view dJones & Puetterlll993h . Weather conditions were poor; 
the thermal background and atmospheric transmission varied 
by 50% throughout the first half of the night. Approximately 
10 individual observations (frames) were obtained at slightly 
overlapping positions using the 12.5 fim filter (12 to 13 /im 
bandpass with > 80% transmission). The chopping secondary 
mirror was driven at 2 Hz with a 30" east-west throw. Each 
frame was observed using the standard mid-IR chop-nod tech- 
nique with two chopping positions plus (on-source) and mi- 
nus (off-source); after chopping with the source in chop-beam 
plus, the telescope was nodded along the chop axis so that the 
object would sit in chop-beam minus, and chopping would 
continue. Each frame was observed for one complete chop- 
nod cycle, yielding a total on-source integration time of 27.6 
s per mosaic frame. The standard stars f3 Peg, /3 And, and /3 
Gem were observed for points-spread function (PSF) determi- 
nation and flux calibration. The images are nearly diffraction 
Umited at 12.5 fim with PSF FWHM = 0.38" a nd a Strehl ra- 
tio of 35%. Data reduction details are given in lShuping et al.l 
dlOOl . The image coordinates are accurate to about 1". 

2.5. Observational results and interpretation 
2.5.1. The CO Outflow Complex 

NRO single-dish observations of accelerated CO emission 
provide constraints on the past evolution of the Cep A out- 
flow complex. While radio jets and HH objects near their 
sources tend to trace the highest velocity and recently ejected 
outflow components with velocities higher than a few hun- 
dred km s"', HH objects located far from their sources and 
shock-excited H2 emission tend to trace moderate speeds of 
10 to about 200 km s"'. CO and other molecules trace the gas 
swept up by secondary interactions with the ambient medium 
and have velocities of a few to tens of km s"'. Thus, CO and 
similar tracers tend to act as calorimeters of the total amount 
of momentum injected into the parent cloud by the faster flow 
components. 

The NRO CO data (Figures [T] |2] and O show a prominent 
outflow complex with complicated structure. In general, the 
Cep A CO emission is collimated along an east-west direc- 
tion at the highest velocities, while at the lowest velocities, 
the CO emission is very poorly collimated and bipolar along a 
northeast-southwest direction. At velocities more than about 
8 km s~' away from the rest velocity of the Cep A core (-11 
km s"' relative to the local standard of rest), the strongest 
outflow component consists of an east-west bipolar flow that 
is blueshifted toward the east (blue component in Figure [U 
and redshifted toward the west. Additionally, an anomalous 
blueshifted lobe appears about 6' west by northwest of HW2 
(upper-right corner of Figure [T]!. A fainter bipolar compo- 
nent, blueshifted toward the northeast and redshifted toward 
the southwest, emerges from the HW2 region at PA. w 35°; 
this is most apparent at moderate velocities (diagonal bipolar 
feature in Figure |2]i. At low velocities of about 2-5 km s"' 
with respect to the centroid velocity, a prominent blueshifted 
lobe emerges from near HW3c (located about 5" south of 
HW2) and extends about 2' east at PA. « 100°. Below, it 
is argued that this feature traces an outflow from HW3c. 

Rather than tracing two independent flows, the brighter 
east-west outflow from HW2 and the secondary, weaker out- 
flow component oriented northeast-southwest at PA. w 35° 
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Figure 3. The 2.12 fim H2 emission (green), and emission from low-velocity CO with — 15 < Vlsr < -13 km s (blue) and —8 < Vlsr < -6 km s (red). 



may trace the walls of a large-scale bipolar cavity that opens 
toward the northeast and southwest. In this interpretation, the 
northeast cavity walls are at P. A. « 100° and 40° and the 
southwest cavity walls are at PA. w 280° and 200°. The ori- 
entation of the IR reflection nebula emerging from the HW2 
region (Figure |5]) provides support for this scenario; its axis 
lies at PA. « 70°, which is aligned with the opening of the 
northeast cavity. 

2.5.2. The Core in the Mid-IR 

Figure |4] shows the Keck 12.5 fim image of the Cep A core. 
The brightest feature is the double-lobed nebulosity surround- 
ing the position of HW2. The two peaks of mid-IR emission 
are separated by about 1.4" along an axis oriented at PA. 
= 45°, which is similar to the orientation of the HW2 radio 
continuum jet. The dark lane that separates the two emis- 
sion peaks has an orientation similar to the circumstellar disk 
surrounding HW2. However, it is unclear whether the dark 
lane is an IR disk shadow produced by a more compact disk 
surrounding HW2, or the gap is actually the disk seen in sil- 
houette. Multicolor imaging is needed to determine the nature 
of the gap. No prominent mid-IR sources are seen at the loca- 
tions of radio sources HW3c and d. However, an extended IR 
nebula is visible about 2" west of HW3c. This feature might 
be illuminated by either HW2, 3c, or 3d. IR polarization mea- 
surements are needed to distinguish these possibilities. 

A faint point source is located at J(2000) = 22:56:18.2, 
H-62:01:44. Although there are no NIR sources at this loca- 
tion, there is an H2O maser spot (H2O-E) located 0.7" from 
the nominal mid-IR source position. 

2.5.3. Radial Velocities of Optical HH Objects 



In contrast to the spectacular 300 km s"' Doppler shifts 
seen toward HH 168 in Cep A West, the long-slit spectra of 
the faint HH objects HH 169 and 174 located east and north- 
east of HW2 show only small Doppler shifts not larger than 
about 60 km s"' in the core of the line profile. All compo- 
nents to the east are blueshifted. The radio proper motions 
of the continuum jet emerging from HW2 indicate veloci- 
ties of order 500 to 300 km s"' within a few seconds of the 
source. Thus, it may seem surprising that shocks located at 
distances of 1-5' away show much lower velocities. The low 
velocities either indicate that the flow east of HW2 is mostly 
along the plane of the sky, or that the ejecta have been de- 
celerated significantly. Future proper motion measurements 
are needed to distinguish these possibilities. However, most 
HH objects exhibit rapid decreases i n their velocities w ith in- 
creasing distance f rom their sources (iReipurth & Ballyll200ll: 
iDevine et al.lll997h . Such deceleration is probably caused by 
the interaction of the ejecta with the ambient medium. The 
structure of the Cep A cloud is such that the outflows emerg- 
ing from the Cep A core propagate relatively freely toward the 
west, but impact dense molecular gas toward the east. Thus, 
it is not surprising that HH 168 has high velocities while the 
eastern shocks associated with HH 169 and 174 have low ve- 
locities. For the rest of the paper, we assume that the space 
velocities of the ejecta east of the Cep A core have a speed of 
about 100 km s"', typical of HH objects located at comparable 
distance from their sources. 

2.5.4. The H2 Outflow Complex 

Figure |5] shows a three-color composite image built from 
the J, H, and Kj broadband filter images. The continuum- 
subtracted narrowband image of the H2 emission in the Cep 
A complex is displayed in Figure |6] and is shown in rela- 
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Figure 4. The Cep A core at 12.5 /im observed with LWS on the Keck 10 m telescope. The locations of the three suspected most massive protostars are shown 
in black circles. HW2 is the unlabeled circle. 



tion to the CO emission in Figures [T}{3] The narrowband im- 
age exhibits three main components: shock emission east of 
the central protostellar cores containing HW2 and HW3c/d, 
bright arcs of shock-excited emission located to the west of 
these cores that is associated with HH 168, and a faint rosette 
of emission to the northwest. The rosette is associated with an 
east-facing globule or small pillar, located about 1.5' north of 
HH 168, which exhibits slightly blueshifted CO emission in 
Figures|2]and|3] The H2 emission is likely to be excited by UV 
radiation from the Cep OB3 association and is thus probably 
fluorescent in nature. The globule contains an mB = 16 mag- 
nitude star that is also called Cep A IRS2. Because there is no 
redshifted emission in this region, the blueshifted CO emis- 
sion is unlikely to be a tracer of an outflow. The CO velocity 
field in the globule may reflect random motions within the 
Cep A cloud, or be the result of UV-photo-heating-induced 
ablation of the globule surface that can accelerate gas to a ve- 
locity of a few km s"' . 

2.5.5. Bipolar Outflow from HW3c: HH 168 and its Counterflow 



The series of bright arcs and H2 bow shocks located 1- 
3' west of HW2 are associated with HH 168. The axes of 
symmetry of most of these shocks indicate a point of ori- 
gin about 10" south of HW2. The sub millimeter and radio 
continuum source HW3c and the bright H2O maser and ra- 
dio source HW3d lie on or near this axis. A dim, 40" di- 
ameter bow shock is located on this axis about 2' east of 
HW3c, directly opposite HH 168. This feature and the visual- 
wavelength part of HH 168 are symmetrically placed about 
HW3c. Thus, HW3c is likely to be the source of these shocks 
and the associated CO outflow co mponents. High-res olution 
cm-wavelength VLA observations (iGarav et al.lll9 96^ reveal 
a chain of radio sources approximately aligned with the axis 
of this flow; these may either trace an ionized jet, or the inter- 
face between a jet and surrounding dense cloud material. 

Figures 2 and 3 show that the bipolar outflow from HW3c 
is associated with low-velocity to intermediate-velocity red- 
and blueshifted gas. The bright jet-like H2 feature located 30" 
east of HW3c is associated with the strongest blueshifted CO 
emission. However, about 1' east, the blueshifted CO lobe 
shifts about 15" south of the flow axis and terminates just 
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Figure 5. Color-combined broadband image of Cep A. J-band emission is shown in blue, H in green, and Ks in red. 
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Figure 6. Continuum-subtracted 2.12 /im H2 emission in the Cep A outflow complex. The successive orientations of a suspected precessing jet from HW2 are 
indicated by red arrows; the position of HW2 is marked with a cross. The oldest (eastward) ejection appears to power HH 174 to the east. The next two ejections 
may be responsible for HH 169. The current orientation of the radio jet emerging from HW2 is at position angle (P.A.) = 45°, indicated by the dashed magenta 
line. Bright HH 168, and the H2 bows at far right, may result from an outflow from HW3c or HW3d (also marked with crosses) along the green axis. A faint H2 
bow marks the opposite lobe of this flow. 



below (south) of the H2 bow shock. Evidently, the HW3c flow 
interacts strongly with ambient cloud material along its south 
wall. However, there appears to be less interaction along the 
north rim of this cavity, perhaps because ambient material in 
this region has been cleared away by the HW2 flow. 

The western lobe of this flow is associated with redshifted 
CO co-located with HH 168. Three bow shocks bright in H2 
lie on the southern side of HH 168; their orientations are con- 
sistent with a flow powered by HW3c. The bow shock mor- 
phologies and low CO velocities indicate that this flow lies 
close to the plane of the sky. It is somewhat unusual that at 
visual and NIR wavelengths, the redshifted outflow lobe is 
brighter than the blueshifted lobe. This feature may be a con- 



sequence of the morphology of the Cep A cloud core that ap- 
pears to be mostly located east of the region containing HW2 
and HW3c. The large range of radial velocities observed in 
HH 168 in visual-wavelength spectra may indicate sideways 
splashing of the high velocity ejecta moving mostly in the 
plane of the sky that has encountered slower-moving mate- 
rial. 

2.5.6. Non-detection of Source W-2 in HH 168 in the IR 

Radio source W is located at the eastern end of the HH 168 
shock complex in Cep A West. The radio emission from 
this object is extended with an elongation similar to that of 
the shocked H2 emission. Based on radio spectral indices. 



8 



Cunningham, Moeckel, & Bally 



iGarav et"aLl d 19961) infer that subcomponent W-2 houses an 
ultracompact internal source, shielded from view in the IR by 
the surrounding molecular cloud, and that the elongated radio 
emission is likely due to a jet launched by this source. Our 
images exhibit co-located, intense NIR H2 emission that co- 
incides with bright optical HH features at the eastern end of 
HH 168, making it unlikely that this region is highly obscured. 
Yet the broadband NIR images do not reveal any point sources 
in this region; thus it seems unlikely that a stellar source ex- 
ists at this location. Rather, radio source W may trace a hard 
shock, with speeds above 400 k ms~'. The detection o f soft 
X-ray emission from this region (iPravdo & Tsuboil2005h pro- 
vides support for this interpretation. Source W is located 
where two outflows, one emerging from HW2 and another 
from HW3c, intersect on the plane of the sky. Thus it is pos- 
sible that these flows are actually colliding. We suggest that 
the radio and X-ray emission from HW may arise from col- 
liding flows emerging from the two most massive protostars 
in Cep A East. This interpretation is generally in line with 
the westward proper motions measured for this radio source 
dRodriguez et al. 2005). 

2.5.7. A Pulsed, Precessing Outflow from HW2 

The blueshifted, eastern lobe of the Cep A outflow contains 
four distinct chains of H2 emitting features, each of which 
terminates in a well-formed bow shock. The four chains ap- 
pear to emerge from the immediate vicinity of HW2. Their 
axes, defined by lines connecting HW2 to the bow shocks at 
the eastern and north-eastern ends of the chains, shift sys- 
tematically clockwise from nearly east-west to northeast- 
southwest. The longest chain, which terminates 4.8' East of 
HW2 at HH 174, has a position angle RA. w 90°. The sec- 
ond chain, which terminates about 3.4' from HW2 at the east- 
ern component of HH 169, has RA. w 80°. The third chain 
terminates at the western component of HH 169, about 2.2' 
northeast of HW2 and has P. A. w 65° . The fourth chain ends 
in a bright but compact H2 bow at RA. « 55° about 1' from 
HW2. The current orientation of the HW2 radio jet contin- 
ues this clockwise migration of outflow orientations and has 
P. A. « 45°. The chains of H2 knots get progressively shorter 
as the axes rotate clockwise toward decreasing PA.. This re- 
markable progression may be an indication that HW2 pow- 
ers a pulsed and precessing jet. Between each major outflow 
ejection event, the jet orientation changes by about 10°-15° 
as seen in projection on the plane of the sky. 

A rough dynamical age for each chain can be estimated by 
dividing the length of each by the velocity of its tip. The ra- 
dial velocities of the HH objects located at the ends of the first 
three chains are low (v < 60 km s"'), indicating that they are 
most likely moving close to the plane of the sky . Fabry-Perot 
imaging of the H2 emission (iHiriart et al.ll2004|) also indicates 
low, but chaotic, radial velocities. On the other hand, the ex- 
citation of the visual wavelength Ha and S II emission of the 
eastern HH objects and comparison with the speeds of other 
HH objects located at similar distances from their sources, 
indicate that shock speeds of order 100 km s"' are not unrea- 
sonable. Future proper motion measurements are needed to 
determine the velocities of the ejecta. Assuming 100 km s"' 
(rea sonable for similar sh ocks exhibiting NIR H2 emission, 
e.g. i Reipurth & Ballvll200 1). the dynamical ages of the four 
eastern H2 chains are 9900, 7000, 4500, and 2100 years, re- 
spectively, indicating that an eruption/ejection event occurs 
approximately every 2500 years. Furthermore, the presence 
of the HW2 radio jet indicates that there is currently an erup- 



tion underway. 

This periodic, geometric progression suggests an underly- 
ing temporal sequence. We propose a model for Cep A East 
in which the accretion disk surrounding HW2, responsible for 
launching a bipolar jet along the disk rotation axis, has period- 
ically changed its orientation. In this picture, the easternmost 
lobe (position angle 90°) constitutes the oldest tracer of past 
jet activity. Since its launch, the accretion disk and associated 
bipolar jet at HW2 have been periodically torqued through a 
series of new orientations. 

Periastron passages of a companion in an eccentric orbit 
that is not coplanar with the disk can readily explain the 
torques required to explain both the impulsive nature of the 
HW2 outflow, and the periodic changes in disk orientation. A 
cluster member in an elliptical orbit around HW2, with a peri- 
astron distance of order the disk radius, will exchange angular 
momentum with the disk. If the binary and disk axes are not 
aligned, the result is a progressive tilting of the disk at each 
periastron passage. The close passage of the orbiting compan- 
ion will also disturb the disk during periastron passage, trig- 
gering an increase in the accretion rate onto the central star, 
and consequent mass loss in the form of a collimated jet. The 
periodic reorientation of the jet axis and disk, and the episodic 
character of the outflow, may be the signatures of a crowded 
and dynamically active environment surrounding HW2. 

The angular extent of the longest shock lobe (position angle 
90°) is 4.8 arcmin, yielding a projec ted length of abou t 1 pc 
given the adopted distance of 725 pc dBlaauw et al.lll959 ). As 
estimated above, the time between disk reorientations (and in 
this model, the orbital period) is sa 2500 years. The mass 
of HW2 is about 15 Mq (Pateletal. 2005, and references 
therein), so a less massive companion must have an orbital 
semimajor axis of w400 AU. The periastron distance must be 
somewhat smaller than this value. As indicated in the Intro- 
duction, Cep A HW2 has at least two companions within the 
required distance, including the "hot core" that may be heated 
by a moderate-mass protostar dMartm-Pinta do et al. 2005D 
and the source of the ex panding maser ring (.Torrelles et all 
l2001btlCurieretani2002l) . 

3. NUMERICAL MODELING 

The interaction of a binary companion and a circumstellar 
disk can cause the disk to p recess (e.g. Papaloizou & TerquemI 
fT99llTerauem et al.l[T999b . This precession may be traced by 
the orientation of the jet, and has been implic ated in the ap- 
pearance of several prec essing outflows (e.g. lEisloffel et alj 
Il996t iDavis et al] 1 1 997h . During repeated encounters of a 
captured companion on an eccentric orbit, t he disk orientation 
move s through several angles impulsively (iMoeckel & Ballyl 
I2OO6I) . rather than smoothly varying as would occur in a 
nearly circular binary. Inspired by the observations presented 
above, we searched the parameter space of encounters simu- 
lated in Mo eckel & B ally (2007a) and selected three combi- 
nations of parameters for further investigation. 

3.1. Method and initial conditions 

We used a modified version of the SPH/A^-body code 
GADGET-2 (Springel 2005) to model interactions between a 
massive protostar and a captured co mpanion. The models an d 
the method are similar to those of Moeckel & Ballvl d2006l) . 
and we summarize them here. A Keplerian disk is set up with 
surface density I](r) cx r~\ and temperature T(r) oc r"'/^. This 
disk is allowed to evolve in isolation until the system has re- 
laxed from its initial conditions. At this time an impactor star 
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Table 2 

Simulation parameters. 



^primary 

Mq 


Mdisk 
Mq 


''disk 

AU 


^impactor 
Mq 


degrees 


^ peri 

AU 


15 


1.5 


350 


5 


30 


150 


15 


1.5 


350 


5 


45 


80 


15 


1.5 


350 


5 


135 


80 



is introduced on a slightly hyperbolic orbit ^ 3 disk radii from 
the primary. The stars are modeled as point masses interact- 
ing with the gas only through gravity, and are assigned an 
accretion radius. Gas that falls into the accretion radius and 
is energetically bound to the star is removed from the calcu- 
lation, and its momentum and mass are added to the star. The 
encounter is integrated until the orbital parameters of the two 
stars have stopped evolving. 

If the orbital timescale of the newly formed binary is much 
larger than the orbital period of the disk, the system is ana- 
lytically advanced to the next interaction. Each gas particle is 
either assigned to the star that it is most bound to or deemed 
ejected from the system, giving us three populations. The 
ejected material is removed from the calculation. In order to 
determine a Keplerian orbit, the gas that remains and the star 
that it is bound to are treated as a point mass at their center 
of mass, and these point masses are analytically advanced to 
the next encounter Upon restarting the hydrodynamic sim- 
ulation, the gas particles are circularized a bout their star us- 
ing a t echnique described in greater detail in lMoeckel & Ballvl 
(|2006|) . The main point of this process is that we avoid a com- 
p utationally expensi v e in tegration to apastron and back. 

'Moeckel & Bally (2007a) performed simulations of en- 
counters between a 20 Mq star with a 2 Mq, 500 AU disk 
and impactors of varying masses, periastron, and inclination 
angle. By inspecting these simulations for disk orientation 
changes of roughly 10°, we selected three impactor parame- 
ters to study the HW2 system, shown in Table |2l We scaled 
the primary and disk masses down to values perhaps more 
suitable for HW2, a 15 Mq primary with 350 AU, 1.5 M© 
disk. The mass of the impactor in all cases is 5 Mq; a 
fairly massive companion is needed to torque the disk through 
~ 10° during a periastron passage. We modeled the disk us- 
ing ~ 1.28 X 10^ particles, and followed each system through 
three encounters. 

3.2. Numerical results 

We consider the orientation of the inner disk Ljisk to be de- 
fined by the summed angular momenta of the gas particles 
within 30 AU of the primary. In order to determine what the 
outflow appearance would be, we construct a crude jet model 
as follows. We assume that the orientation of the outflow is 
the same as the inner disk orientation, Vjet = i-disk- By assum- 
ing a constant- velocity jet launched in this direction, we can 
construct a model of what the jet from this interacting system 
will look like. Of the three cases considered, the two pro- 
grade simulations (/ = 30° and / = 45°) produced jets that are 
similar in appearance to the observed system. The retrograde 
encounter, / = 135°, showed a shift in orientation similar to the 
other two for the first passage, but the subsequent passages did 
not resemble the actual jet. We believe this is in part due to 
the rapid disk destruction associated with retrograde passages 
having small rperi/raisk ( Moeckel & Bally 2006). 

Plots of the jets through the first three passages from the 
two prograde simulations are shown in Figure |7] They have 



been rotated about the initial jet axis in order to most closely 
match the observed appearance of the HW2 system, and the 
initial jet axis lies on the x-axis to match the orientation of the 
jet lobes in Figure |6] The simulations yield the elapsed time 
of the encounters, and at an assumed distance of 725 pc, the 
distance to the oldest bowshock is approximately 1 pc. By 
choosing the velocity of the jet to match this distance, we find 
jet velocities of 88 km s"' for the 30° case, and 139 km s"' 
for the 45° case. The boxes on the figure indicate where bow- 
shocks would be expected, and are labeled with the age of the 
flow at that point. While both simulations capture the qual- 
itative orientation of the jet, the case with / = 45° appears to 
match the observed jet orientation slightly better. 

We stress that we are not claiming that the parameters stud- 
ied here reflect the reality of the HW2 region's dynamics. The 
uncertainties in the masses, radii, and even number of objects 
impede attempts to model the actual system, and the range 
of parameters that produce similar outflows is large. Rather, 
these simulations should be viewed as a proof-of-concept ex- 
periment. Using plausible values for the masses of the com- 
ponents, we can generate outflow orientations similar to those 
observed, via a series of postcapture interactions between a 
massive protostellar system and a moderate-mass companion. 

4. DISCUSSION 

The two most massive protostellar objects in Cep A, HW2 
and HW3c, each power separate outflows. HW2 appears to 
have been responsible for at least four quasi-periodic ejec- 
tions that produced collimated flows visible in NIR H2 emis- 
sion in the eastern lobe of the Cep A outflow complex. While 
the first eruption launched a flow nearly due east, subsequent 
ejections indicate that the source orientation rotated north by 
about 10° between events. The current orientation of the radio 
jet, rotated clockwise by 45° compared to the ejection axis of 
the first event, continues this trend. Radial velocity measure- 
ments of the visual wavelength emission associated with HH 
169 and 174 located at the tips of three of the outflow axes in- 
dicate low radial velocities and dispersions, an indication that 
the outflow axis probably lies close to the plane of the sky. 

The bright NIR reflection nebula emerging from the vicinity 
of HW2 has an axis close to the last two ejection events. Mil- 
limeter wavelength CO emission toward the eastern lobe of 
the Cep A outflow complex reveals two ridges of blueshifted 
emission at P.A. ^ 100° and 35°, approximately symmetri- 
cally placed about the axis of the reflection nebula. A pair of 
redshifted ridges of CO emission are located on the opposite 
side of HW2 toward the southwest. These CO features may 
define the walls of an outflow cavity created by the outflow 
over time. The CO outflow from HW3c is superimposed on 
the east-west rim of this suspected cavity. 

The data presented above are interpreted in terms of a 
pulsed, precessing jet emerging from Cep A HW2. A numer- 
ical model is used to show that a pulsed, precessing jet can be 
the result of disk perturbations produced by a companion star 
in an eccentric orbit that is not co-planar with the disk. This 
scenario is the most li kely outcome o f the capture-formed bi- 
nary proposed by (Moeckel & Ballv 2007a b). We note that 
[Narayanan & Walker (1996) also suggested periodic outflow 
from Cep A, although the dynamical timescale they deter- 
mined for the launch of successive CO shells is ^ 1 .6 x 10^ yr, 
two orders of magnitude greater than the ^ 2500 yr timescale 
suggested by the higher velocity shocked H2 gas presented 
here. 

The second outflow from the Cep A core region appears to 
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Figure 7. Locations of jet ejecta from two simulations, each as viewed at a single time. The solid lines mark projected distance from outflow source, with 
older material at far left and recent ejecta at far right. Dotted lines indicate successive jet orientations; in the absence of disk precession, the jet would remain 
horizontal. The jet velocities have been scaled to give a 1 pc distance to the oldest shock, and the initial jet orientation has been chosen for comparison to Figure 
[6] Boxes indicate where shocks at the ends of the outflow lobes would be expected, and are labeled with the age of the flow at that point. 



emerge from slightly south of HW2. This flow is associated 
with a low-velocity blueshifted lobe of CO emission that ter- 
minates about 1 .5 to 2' to the east in a faint but large H2 bow 
shock. HW3d is associated with strong H2O masers that are 
located at the base of a chain of radio sources at position angle 
P.A. « 105° dTorrelles et alll 99811200 fab fliat has the spectral 
index and morphology of a thermal radio jet that has the same 
orientation as the suspected large-scale outflow from HW3c. 

A third young star located southwest of HW2 may also 
power an outflow that produced the radio sources HWlb, 
3a, and 3b. HW3b is the third source in Cep A to con- 
tain masers (iTorrelles et al.lll998ll2001a|). The sub millimeter 
source SMA4 (Broganet al. 2007) is located on the axis of 
this chain close to HW3b. The chain of radio sources has a 
position angle P.A. « 60-65° (or 240 to 245°), close to the 
orientation of the H2 jet associated with the third suspected 
eruption of HW2. Thus, an alternative to our model is that 
this collimated component east of the Cep A core is a steady 
jet associated with a protostar located in the HWlb, 3a, and 3b 
chain. If SMA4 is the source, then the chain of radio sources 
is the counter-jet associated with the prominent H2 feature. 
However, this model does not explain the orderly progression 
of terminal bow-shock orientations and distances discussed 
above. Furthermore, the chain of radio sources HWlb, 3a, 
and 3b lies in the "line-of-fire" of the current HW2 jet. Thus, 
it is possible to interpret these features as shocks where the 
HW2 jet rams part of the Cep A molecular cloud, or the out- 
flow from HW3c. 

No NIR continuum source is detected at the location of ra- 
dio source W at the eastern end of HH 168. The presence 
of radio continuum and X-ray emission from this portion of 
HH 168 may indicate that this region is being impacted by 
shocks with speeds of at least 400 km s"' . We propose that 



HH 168 marks the location where the HW3c outflow collides 
with the outflow from HW2. The location of HH 168 indi- 
cates that the currently shocked gas is mostly excited by the 
debris from the third and fourth eruptions of HW2 with the 
flow form HW3c. Ejecta from the fifth eruption, which is 
traced by the HW2 radio continuum jet, has not yet reached 
the interaction region. If HW3b also drives its own jet with 
an orientation similar to the third eruption from HW2, it may 
also contribute to the excitation of HH 168. 

The complex shock morphology may be the result of tur- 
bulent mixing of the colliding flows. The north-to-south den- 
sity gradient in the ambient molecular cloud may contribute 
to the bending of this shock complex toward the north. The 
CO emission, which is both red- and blueshifted at the loca- 
tion of HH 168, may also be deflected along our line of sight 
by this gradient to produce the blueshifted CO lobe farther 
downstream. 

Collisions between outflows and clouds have been in- 
voked to explain the complex structure of some Herbig- 
Haro objects such as HH 110 in Orion where the HH 270 
jet sl ams into a small m olecular cloud in the L1617 com- 
plex (iLopez et alJ 120051) . An embedded IR source in this 
cloud , IRAS 05487-H0255, dri v es a jet a nd molecular ovX- 
flow jReipurdi & Olberd [19911: iLee et all I2OOOD . Kajdic et 
al. (2009, in preparation) propose that HH 110 is formed 
by the collision of the HH 270 flow with the flow from 
IRAS 05487H-0255. Th eoretical models of je t -cloud interac- 
tions were presented bvlCanto & Ragal (ll996h:lRaga & Cantol 
(ll996l) : ldeGouveiaDal Pinol (ll999^:'Raga et al.' ('2002'). Jet- 
jet collisions were modeled by Cunningham et al. (2006). In 
Cep A, HH 168 may be excited by the collision between a 
pair of outflows and be additionally deflected by the density 
gradient in the cloud. Observations show that the sources of 



Precessing Jet in Cep A 



11 



the HW2 and 3c flows are physically close. Furthermore, both 
sources appear to drive flows that are close to the plane of the 
sky. These factors make a collision of flows reasonably prob- 
able. The absence of obvious shocks southwest of the Cep 
A core supports the view that ejecta from HW2 are deflected. 
Still, the presence of low-velocity CO emission southwest of 
HW2 and HW3c indicates that this part of the cloud has been 
impacted by outflow activity. A possible explanation for this 
dichotomy is that, in the distant past, only the outflow from 
HW2 was active, and that the deflection of this flow by the 
flow from HW3c is a relatively recent phenomenon. 

5. CONCLUSIONS 

We present new narrowband NIR images of the the Cep A 
outflow complex in the 2. 12 /im S(l) line of H2. These images 
provide evidence for two major outflows from the two most 
luminous members of the Cep A outflow complex. 

The radio source HW3c appears to drive a collimated east- 
west outflow with a redshifted lobe emerging toward P.A. « 
275° that may be responsible for much of the emission associ- 
ated with HH 168 in Cep A West. The blueshifted lobe of this 
outflow appears to terminate in a large but faint H2 bow shock 
located about 1 .5-2' east of HW3c. The H2 bow and the the 
brightest portion of HH 168 are placed symmetrically about 
HW3c. The three H2 bow shocks located along the south- 
ern portion of HH 168 have orientations consistent with being 
driven by this source. 

The radio source HW2, the most luminous and massive pro- 
tostar in the Cep A complex, appears to drive a precessing and 
pulsed jet. The H2 images reveal four distinct jet axes. The 
longest and oldest outflow lobe terminates in HH 174 about 
5' due east of HW2. The second and third lobes are pro- 
gressively shorter, have rotated clockwise by about 10-15° 
between ejections, and terminate in the western and eastern 
components of HH 169. The fourth ejection is the shortest 
lobe and was ejected toward the northeast, but has no visual 
wavelength HH objects associated with it, probably due to the 
high obscuration. The radio continuum jet may trace a cur- 
rent period of activity and continues the trend of clockwise 
rotation of the outflow axes. Thus, over the past 10"* years, 
the HW2 outflow appears to have precessed by nearly 45° in 
a clockwise direction as seen on the plane of the sky. Assum- 
ing an average flow speed of 100 km s"^ HW2 undergoes an 
eruption about every 2500 years. 

Most massive stars are born in clusters, and Cep A is no 
exception. In an environment with a high stellar density, in- 
teractions between a massive star and its siblings become rel- 



atively common. Cep A HW2 is surrounded by a disk at 
least several hundred AU in radius and with a mass of order 
1 Mq. This disk can serve as a dissipative environment, and 
any star passing within about a disk radius has a high prob- 
ability of being captured into a highly eccentric orbit. The 
more massive the intruder, th e higher its probability of cap- 
ture (iMoeckel & Ballvll2007allbl) . This mechanism provides a 
natural explanation for the presence of a moderate-mass com- 
panion, indicated in the HW2 s ystem by a hot core 0.6" (40 
AU in projection) east of HW2 dMartm-Pintad o et al. 2005). 

It is proposed that Cep A HW2 has a moderate-mass com- 
panion in an eccentric orbit whose orbital plane is inclined 
with respect to the circumstellar disk surrounding HW2. Pe- 
riastron passages of the companion may perturb the disk, 
and drive accretion onto the central star, and produce quasi- 
periodic episodes of collimated mass ejection. The non- 
coplanar eccentric orbit also applies a torque to the disk, 
changing its orientation. We investigate this hypothesis us- 
ing an SPH code and demonstrate that such interactions can 
result in the type of disk orientation change proposed to have 
occurred in Cep A. Modeling indicates that a companion in 
a prograde orbit inclined with respect to the HW2 disk by 
about 40° fits the observations best. Because an impactor 
mass greater than the disk mass is needed to torque the disk 
through an appreciable angle, the presence of one or more 
low-mass sources in the vicinity do not affect this model. 

The bright shock complex HH 168 may be excited by the 
collision of two or more outflows emerging from the Cep A 
cloud core located about 90" to the east. The western end of 
HH 168 may trace the chaotic mixing resulting from the col- 
lision of the steady (in orientation) flow from HW3c and the 
third eruption of HW2; the eastern end of HH 168 may trace 
the collision of the flow from HW3c with the fourth eruption 
of HW2. If SMA4 also drives a flow, it may also participate 
in this interaction. 
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